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Stochastic Beta Relaxation (SBR) is a model for the dynamics of glass- forming liquids close to 
the glass transition singularity of the idealized mode- coupling theory (MCT) that has been derived 
from generic MCT-like theories by applying dynamical field-theory techniques. SBR displays a rich 
phenomenology common to most super-cooled liquids. In its simplest version it naturally explains 
t-wo prominent features of the dynamical crossover: the change from a power-law to exponential 
increase in the structural relaxation time and the violation of the Stokes-Einstein relation between 
diffusion and viscosity. The solution of the model in three dimensions unveils a qualitative change 
at the crossover in the structure of dynamical fluctuations from a regime characterized by power-law 
increases of their amplitude and size to a regime dominated by strong Dynamical Heterogeneities: 
rare regions where dynamics is relatively much faster than in the rest of the system. While the 
relaxation time changes by orders of magnitude, the size of these regions does not change significantly 
and actually decreases below the crossover temperature. SBR cannot sustain too large fluctuations 
and could fail well below the crossover temperature. There it could be replaced by non-conventional 
activated dynamics characterized by elementary events with intrinsic time and length scales of an 
unusual large (but not necessarily increasing) size (mesoscopic vs. microscopic). 

PACS numbers: 64.70.Q 


I. INTRODUCTION 

Mode-coupling-Theory provides a good qualitative and 
quantitative description of the initial dynamical slowing 
down of super-cooled liquids but its main problem is that 
it predicts dynamical arrest at a temperature Tc where 
a crossover is observed instead [T]. Recently it has been 
shown that the solution to this problem may come by 
treating the singularity as a genuine phase transition by 
means of perturbative field-theoretical methods [2] • This 
is surprising because in general perturbative methods are 
unable to remove a singularity, however one can show 
that in the case of MCT the perturbative loop correc¬ 
tions are the same of those of some dynamical stochastic 
equations called stochastic-beta-relaxation (SBR) equa¬ 
tions in [2]. If studied perturbatively both MCT-like 
theories and SBR display dynamical arrest at all orders, 
but the SBR equations can be also solved explicitly (be. 
non perturbatively) showing that the transition is instead 
changed into a crossover due to non-perturbative effects 
that can be clearly identified. 

The result is rather intuitive: on the time-scale of the 
/3-regime (where by definition the density-density corre¬ 
lator remains close to a plateau) dynamics according to 
the SBR equations is described by the very same MCT 
critical equations, the only difference being that the tem¬ 
perature fluctuates randomly between different regions of 
the system. As a consequence, even if the global temper¬ 
ature is near Tc or below, there are regions of the system 
in the liquid phase and they destabilize the glass phase 


predicted by ideal MCT and restore ergodicity. 

In a previous publication [3] we have studied a sort of 
schematic version of the model where different regions are 
uncorrelated and behave independently. The study of the 
model requires elementary computations but, notwith¬ 
standing its simplicity, displays many features typically 
observed in super-cooled liquids and allows to understand 
them in an intuitive way. In particular when supple¬ 
mented with the assumption of time-temperature super¬ 
position this simplified SBR allows to obtain predictions 
on the a regime that are in remarkable agreement with 
the known phenomenology of various quantities, includ¬ 
ing the a-relaxation time, the Diffusion constant and the 
thermal susceptibility. The main limit of the simplified 
model is that it does not allow to study length scales 
which requires the solution of the full SBR equation in 
finite dimension. 

In this paper we will discuss the solution of the SBR 
equations in three dimensions. From the solution one 
can draw a rather comprehensive description of the glass 
crossover that we will sketch in the following. At any 
point in space we associate some local field that quanti¬ 
fies somehow the mobility i. e. the time rate with which 
that portion of the system decorrelates from the initial 
configuration. The amplitude and the size of the fluctua¬ 
tions of the local mobility allows to quantify and discuss 
the notion of Dynamical Heterogeneities. The solution 
of the SBR equations suggest that 

• approaching Tc from above the mobility field (as¬ 
sociated to the function B{x) in the following) de- 


2 


creases in average value (it would be zero in ideal 
MCT at Tc) while its local fluctuations increases 
both in size and in amplitude. More precisely the 
process has the features of a second-order scale- 
invariant phase transition and in particular it is 
characterized by an increasing dynamical correla¬ 
tion length. However the increase of the relative 
fluctuations is not very pronounced approaching 
from above and it seems not appropriate to talk of 
Dynamical Heterogeneities in the sense used by ex¬ 
perimentalist. 

• Close to Tc, there is a change and two important 
things happens: i) dynamical arrest is avoided and 
ii) the relaxation time starts to grow much faster, 
from power-law to exponential-like. 

• The structure of dynamical fluctuations also dis¬ 
plays a qualitative change below Tf.'. overall the 
dynamics continue to slow down dramatically but 
it is dominated by rare regions that are relatively 
much faster than the rest of the system. Dynamical 
Heterogeneities are the hallmark of the dynamics in 
this regime. 

• Below Tc dynamics is slow not because the faster 
regions are large, but rather because they are rare 
and therefore the amplitude of fluctuations in¬ 
creases. Actually the size of these regions shrinks 
while decreasing the temperature below Tc and the 
dynamical correlation length decreases. In general 
the change in the correlation length is not dramatic 
and decouples from the increase of the relaxation 
time. The key point is that the structure of the fluc¬ 
tuations changes from scale-invariant-like above Tc 
to activated-like below Tc. Correspondingly quanti¬ 
ties that scales similarly above Tc decouples leading 
e.g. to deviations from the Stokes-Einstein rela¬ 
tionship (SER). 

In the next section we will discuss how the above pic¬ 
ture emerges from the numerical solution of the SBR 
equations in 3D and in the final section we will discuss 
the results. 


II. A THEORY OF THE GLASS CROSSOVER 
A. Stochastic Beta Relaxation 

MCT and similar mean-field theories lead to the pre¬ 
diction that near the critical temperature the density- 
density correlator has the following behavior on the time 
scale of the /3-regime Tp: 

=F(fc) + G(t)H(fc) (1) 


function G{t) obeys the well-known MCT equation for 
the critical correlator: 

d f* 

a = -\G'^{t) + — / G{t - s)G{s)ds 
dt Jo 

where the separation parameter a is negative at high tem¬ 
peratures (low pressures) and vanishes at the MCT sin¬ 
gularity: 

a(xTc — T,a(xp — Pc 

The solution of the above equation is such that G{t) goes 
to minus infinity at large times in the liquid phase accord¬ 
ing to the so-called Von Schweidler’s law: 

lim G(t) = for cr < 0 (liquid) 

i—>-oo 

and it goes instead to a constant in the glassy phase, 
signalling that the correlator remains blocked near the 
ergodicity breaking parameter F'^{k): 

lim G(t) = -y/tT/(l — A) for (T > 0 (glass) 

t—>-oo 

The exponent b is expressed in term of the exponent pa¬ 
rameter A through: 

^ r^(i-a) ^ r^(i + 6 ) . , 

r(l-2a) r(l-H26) 

Where a is the exponent controlling the small time be¬ 
havior of G{t) ~ l/t“both in the liquid and glassy phase. 
The critical equation is valid provided that G{t) is small 
and this condition defines the j3 regime. When G{t) be¬ 
comes 0 ( 1 ) we enter the a regime, whose time-scale can 
therefore be obtained as: 


Ta oc B{a) (3) 

The prefactor B{a) of the term vanishes in ideal 
MCT approaching cr = 0 as |tT|'>'^ (where 7 = l/(2a) -I- 
1/(26)) and as a consequence diverges as Tq ~ |cr|“'>' . 

SBR can be viewed as an extension of the MCT equa¬ 
tion for the critical correlator with random fluctuations of 
the separation parameter. According to it in the (3 regime 
equation Q continues to hold and only the epxression of 
the critical correlator G{t) is different. One must consider 
a field g{x, t) that is a local version of the correlator and 
that obeys the following equation: 


d 

a+s{x) = —aV'^ g(x,t)—Xg'^(x,t)+— / g{x,t—s)g(x, s)ds 

at Jo 

(4) 

where the field s{x) is a time-independent random fluc¬ 
tuation of the separation parameter, Gaussian and delta- 
correlated in space: 


[s(a:)]=0, [s{x)s{y)] = Aa'^ 6{x - y) (5) 


The total correlator is obtained as the integral over space 
averaged over the random fluctuations: 


where the bold character accounts for the case of mix¬ 
tures of particles in which the correlator is a matrix. The 


Git) = [^ / 9ix,t)da 


( 6 ) 
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FIG. 1. Top: Pictorial representation of the solution of SBR 
equations in finite dimension for a given realization of the ran¬ 
dom s{x). (generated from actual solution of a ID system). 
At large times the solution converges to the form 
thus inducing a mapping between the realization of the ran¬ 
dom s(x) and a positive function B(x), the B-profile. Bottom: 
Plot of g{x,t)/t^ vs. X for increasing times t: increasing the 
time the curves converge to the B-profile. 


Thus SBR introduces two novel parameters in the de¬ 
scription of the (3 regime: the variance Acr^ of the ran¬ 
dom fluctuations and the length-scale of spatial corre¬ 
lation which are controlled by the coefficient a in front 
of the nabla term. The fact that SBR can be derived 
from MCT-like theories at criticality implies that quanti¬ 
tatively these parameters can be computed within MCT 
similarly to A and a [2], for instance the parameter a 
can be extracted from Inhomogeneous MCT [4]. Differ¬ 
ent estimates can be obtained in native field-theoretical 
approaches [S]. 

B. The B-profile 

In the simplified model the dynamics in each region de¬ 
pends solely on the given value of the local temperature 
and this allows to understand two important features of 
SBR [^1311^. First of all one can see that at any temper¬ 
ature, even deeply below Tc there will be liquid regions 
because of fluctuations of s(a;), therefore dynamical ar¬ 
rest is avoided. Second, the typical region is liquid above 
Tc but is frozen below and correspondingly the fluctu¬ 
ations s{x) ^ —a required to have a liquid region become 
increasingly rare below Tc; this determine a crossover in 
the growth of the relaxation time from power-law to ex¬ 
ponential. We have studied the SBR equation in 3D 
numerically solving for g{x,t) at given s{x) and veri- 


FIG. 2. Normalized B-profile B{x)/B on a plane sliced from 
a cubic box for different values of a. 


fied that the above description remains essentially valid. 
In particular both for negative and positive values of a 
the solution g{x,t) escapes to minus infinity {i.e. leaves 
the plateau) at all points x. This corrects a disturbing 
pathology of the simplified SBR equations where there 
are always regions that remains blocked near the plateau 
and never decay. More precisely g{x,t) exits from the 
plateau with the very same Von-Schweidler law with a 
space dependent constant B{x), (see hg. [^) 

g{x,t) ~ —Bcr{x)t^ for all a 

Therefore at any value of a the SBR equations induce a 
non-trivial mapping between the realization of the fluctu¬ 
ations s{x) and a positive function B{x) (called B-profile 
in the following). In practice we extract the B-profile by 
solving the equations up to times large enough to be in 
the asymptotic regime where g{x,t) « —B{x)t^. The 
presence of the gradient term in the full SBR equation 
leads to the disappearance of a clear distinction between 
liquid and glassy regions as all regions become liquid, 
nevertheless when we switch on the gradient the regions 
that were glassy will be characterized by a very small 
value of B{x) compared to the rare liquid regions where 
B(x) is relatively much larger. The B-profile will be our 
main focus in the following, indeed it allows a discussion 
in a clear and compact way of many dynamical quanti¬ 
ties including the a-relaxation time, diffusion coefficient, 
correlation length and Dynamical Heterogeneities. 

It is illuminating to directly inspect the B-profile above 
and below the critical temperature. In figure ([^ we plot 
the normalized B-profile B{x)/B on a plane sliced from 
a cubic box. The peaks(valleys) in the B-prohle corre¬ 
spond to regions that are decorrelating from the initial 
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condition faster(slower) than the average and where lo¬ 
cal dynamics is also faster (slower). The height and size of 
the peaks allows to characterize dynamical fluctuations. 
We see that for T > Tc { i.e. cr < 0) fluctuations are not 
very pronounced and the system appear homogeneous, 
instead below Tc more and more peaks disappear and as 
a consequence the few that are left tend to be much more 
pronounced. Overall the average value B decrease mono- 
tonically with decreasing the temperature (it changes by 
orders of magnitude for the values of a of the plot) but 
the dynamics becomes instead dominated by rare regions 
that are relatively much faster than the typical region. 
Thus the evolution of the i?-profile encodes a qualitative 
change in the structure of dynamical fluctuations marked 
by the appearance of Dynamical Heterogeneities at the 
crossover temperature. On the other hand one can im¬ 
mediately see that the size of the peaks does not change 
significantly above and below Tc. In the following sec¬ 
tions we will analyze the H-profile more carefully and we 
will see in particular that the correlation length is actu¬ 
ally slightly decreasing below Tc- 

C. Viscosity and DifFusivity 

The viscosity r] and the diffusivity D can be associated 
to different averages of the H-profile. Following the same 
matching arguments valid in MOT we can associate the 
a time scale to the coefficient of in the total correlator 
G'(t), this leads naturally to: 

f] = Ta ^ J B{x)dx^ (7) 

On the other hand, as discussed in [3], within SBR it 
is natural to consider a local relaxation time Ta(x) ^ 
B(x)~^''^ defined as the time scale where —B{x)t^ is 
0(1). One would therefore associate the diffusivity to 
the inverse of this local relaxation time: 

D^^jB^^\x)dx (8) 

In figure ([^ we plot and I? as a function of cr to¬ 
gether with exemplary experimental data, obtained by 
Schneider et al. El H] from dielectric spectroscopy on 
propylene carbonate. The behavior of this quantities is 
similar to what has been found in the simplified model. 
Let us make a few comments on the figure: 

• the divergence of Tq, at cr = 0 of ideal MCT is 
avoided but the growth rate changes from power- 
law to a more pronounced exponential-like growth: 
Tc is avoided but marks a crossover. 

• The above property should be assessed remember¬ 
ing that there is no ad hoc assumption on activa¬ 
tion processes in the derivation of SBR, rather the 
initial assumption is i.e. that Tc marks a genuine 
phase transition. 



- 0.04 - 0.03 - 0.02 - 0.01 0.00 0.01 0.02 

cr 


FIG. 3. Top: Ta vs. a from SBR in 3D. Inset: D vs. o 
in 3D. Symbols are dielectric-spectroscopy data for propylene 
carbonate from Lunkenheimer et al (2000). Dashed lines are 
MCT asymptotes. Bottom: The correlation length (de¬ 
fined as half-width at half-maximum of r(r)) as a function of 
a. 



-0.04 -0.03 -0.02 -0.01 0.00 O.OI 

FIG. 4. Top: Ta vs. cr from SBR in 3D. A power-law fit 
of the SBR expression would suggest a critical temperature 
dehnitively smaller than the actual value cr = 0. Quantita¬ 
tively the value of the shift depends on the parameters Act 
and a. 


• The comparison with the data show that the de¬ 
scription provided by SBR could work in an ex¬ 
tended range of temperature, although it depends 
on a few number of parameters. 

• Used as fit functions the SBR expressions for D 
and rj could help reconcile different estimates of Tc 
based on asymptotic behavior. 
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• SBR leads to an apparent shift of Tc to lower val¬ 
ues when fitted with ideal MCT power laws and 
this could provide an explanation of the well-known 
shift of the quantitative prediction of MCT with 
respect to simulation data [T], see fig. Q. 

• In the figure the SBR expression for the was 
used essentially as a fit function to the data but 
the actual computation of the parameters Act and 
a is feasible for many systems that can be simulated 
numerically. These computations are left for future 
work. 

As a technical remark we note that at given value of 
A the SBR equations depend in principle on three pa¬ 
rameters a, CT and Act. However it suffices to solve it for 
fixed values of two of them. The natural way to do that 
is to consider fixed values of Act and a and vary a. The 
general solution can then be obtained from the reference 
solution by means of appropriate rescalings. The numer¬ 
ical data shown here correspond to the following choice 
of the parameters: 

a = .2, Act = .1 , A = .75 , (9) 


D. Non-monotonous Correlation Length 


The spatial fluctuations of the dynamics in the late /3 
regime are conveniently encoded by the spatial fluctua¬ 
tions of the R-profile. We introduce the self-correlation 
of the H-profile as: 

r(r) = :^ J B{x)B{x + r) dx (10) 


The numerical solution shows that while the absolute 
value of r(r) changes by orders of magnitude upon cross¬ 
ing CT = 0 (it is related to Tq,), its shape and length-scale 
vary much less. In order to focus solely on the space 
dependence we introduce the normalized self-correlation: 


rn(?') 


r(r) -r(oo) 
r(0) - r(oo) 


( 11 ) 


We observe that in 3D r„(r) has a bell-shaped form with 
rapidly decaying tails both below and above ct = 0. Quite 
interestingly it turns out that the width of r„(r) has a 
non-monotonous behavior with ct: in the bottom of figure 
([^ we plot the half-width at half-maximum of r„(r) as a 
function of ct and use it as a definition of the dynamical 
correlation length. The figure demonstrates what we an¬ 
ticipated in the introduction: the dynamical correlation 
length increases approaching ct = 0 from the liquid phase, 
it saturates to a maximum at ct = CT^ax slightly below 
Tc and then it decreases again. In the same region the 
relaxation time increases instead by orders of magnitude 
implying that near ct = 0 the correlation length decouples 
from the relaxation time. For instance considering the 
interval ct = [—.01, .015] the correlation length increases 


and decreases again with an approximate 25% excursion 
while the relaxation time increases by almost 8 orders of 
magnitude. 

It is tempting to put these results in connection 
with recent observations of non-monotonous correlation 
lengths in numerical simulations [9] and experiments 
m- These observations however should be contrasted 
with recent measurements mm of dynamical corre¬ 
lation lengths for the same systems obtained with the 
more standard methods of Refs. [IMl]. In this case 
a monotonous behavior was instead observed, albeit 
displaying evidences of saturation towards a maximum 
mm- With regard to this open issue we note that, as 
discussed in [5], SBR is valid as it is near the crossover 
region. In particular the assumption that the separation 
parameter ct is the sole temperature-dependent quantity 
(a typical assumption for a genuine phase transition) 
could be too strong well below T^. On the other hand ac¬ 
cording to fig. ([^ the correlation length does not change 
too much with the temperature, and one cannot exclude 
the possibility that the temperature dependence of the 
parameters a and Act alters the non-monotonous behav¬ 
ior in actual systems. Qualitatively however the scenario 
of fig. ([^ would remain the same: the amplitude of fluc¬ 
tuations (the height of the peaks) increase considerably 
while their size does not change significantly and this 
should be considered the essential feature of SBR inde¬ 
pendently on weather this size is actually increasing or 
decreasing. 


E. Dynamical Heterogeneitites 

Given that the correlation length displays a symmet¬ 
rically decrease far away from Tc both above and below 
one may ask what is actually happening at the crossover. 
As we saw before the answer is that there is a dramatic 
change in the structure of dynamical fluctuations above 
and below the crossover temperature as can be seen by 
direct inspection of the H-profiles. Let us now examine 
figure thoroughly. We can distinguish two regimes 
above and below the crossover temperature. Above the 
critical temperature (ct < 0 ) the normalized R-profile has 
fluctuations of order 0(1) around its average value which 
is one by definition. The profile is characterized by peaks 
whose width is of the order of magnitude of the correla¬ 
tion length. Increasing ct towards ct = 0 the profile does 
not change much, although we observe slightly larger fluc¬ 
tuations in the height of the peaks. In this regime the 
profile evolves much as in a second-order phase transition: 
the correlation length increases and so do the amplitude 
fluctuations. Here the correlation length carries relevant 
information: the system is essentially scale-invariant in 
the sense that the profiles at different values of ct < 0 look 
the same once space is rescaled proportionally to the cor¬ 
relation length. However at higher values of ct (ct = .006) 
the R-profile starts to change qualitatively: the height of 
the typical region decreases as more and more of the low- 
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FIG. 5. Viscosity-Diffusivity parametric plot from SBR in 3D. 
Dashed line indicates D oc ri~^ (Stokes-Einstein relation), a 
large circle marks T = T^. Smaller symbols: experimental 
data for o-terphenyl from Ref. |18| (circles: tracer diffusion at 
T > Tc’, squares: diffusion of flourescent ACR dye; diamonds: 
TTI dye); simulation results for a harmonic-sphere mixture 
(triangles, Ref. |12|i. 


est peaks disappear, conversely the few peaks left increase 
their relative height because they carry all the weight. As 
a result at a = .012 we have a completely different land¬ 
scape characterized by rare regions where the dynamics 
is relatively much faster {B{x)/B ^ 1 at the peaks) then 
the surrounding slowly-moving regions {B(x)/B ^ 1 in 
the typical region). Quantitatively we see that while 
fluctuations are 0(1) at cr = —.006, for a = 0.012 the 
normalized i3-profile is 2 — 3 orders of magnitude larger 
in the rare fast regions with respect to the neighboring 
slow regions. From figure ([^ we also see that the size and 
shape of the peaks does not change much, which result in 
the fact that the correlation length does not change sig¬ 
nificantly for the values of a considered (as seen in figure 
(§)• 

These features suggest that we should consider care¬ 
fully the role of the correlation length at the glass 
crossover and the very same notion of Dynamical Hetero¬ 
geneity (DH). Although dynamical fluctuations grow in 
both regimes they have power-law increase and remains 
relatively small in the first regime while they grow signif¬ 
icantly in the second regime where they have essentially 
an activated nature. Therefore only in the second regime 
we should actually talk of Dynamical Heterogeneities as 
defined from experiments m and as a consequence we 
should conclude that within SBR DH are not intrinsically 
associated with an inereasing correlation length. 

It is to be expected that the qualitative change in 
the structure of the dynamical fluctuations from scale- 
invariant-like to activated-like can be detected by com¬ 
paring observables that are associated to different aver¬ 
ages oi B{x). In figure (§ we plot parametrically the 
viscosity and the diffusivity computed according to ex¬ 


pressions Q and They obey the Stokes-Einstein 
relationship (SER) Drj oc I for temperatures above 
while below the crossover (the black dot) there is a vio¬ 
lation of the SER. Once again this is a typical feature of 
glassy systems as the exemplary data demonstrates. 


III. DISCUSSION 

In order to assess the properties of SBR described be¬ 
fore one should bare in mind that it is not a phenomeno¬ 
logical theory and there is instead a non-trivial math¬ 
ematical connection between MCT and similar micro¬ 
scopic theories characterized by an ideal glass transition. 
As we said in the introduction this connection is rather 
unexpected because it is based on the application of per¬ 
turbative methods to an effective dynamical field theory 
which coincides at the tree level with the equation for the 
critical correlator of ideal MCT. 

The dynamical field theory from which SBR was de¬ 
rived is closely related to a statie replicated field theory 
that was associated to the glass problem long ago . 

In particular a static stochastic equation was first derived 
from the replicated field theory in [^. One should be 
aware that this is the same (static) field theory that lies at 
the heart of the Random-Fisrt-Order-Transition (RFOT) 
theory [23] . As it is well known, RFOT claims to include 
the early stage of vitrification described by MCT phe¬ 
nomenology but is definitively more focused on the exis¬ 
tence of an ideal glass transition below the calorimetric 
glass transition Tg and put emphasis on the fact that 
the corresponding dynamical slowing down is accompa¬ 
nied by a diverging correlation length. SBR is instead a 
theory of the glass crossover and it is not clear a priori 
how deep in the super-cooled regime it provides a good 
description. Thus the two theories are not necessarily in¬ 
compatible because they describe different temperature 
regimes. On the other hand clearly if one were to specu¬ 
late on the deeply supercooled regime starting solely from 
the SBR description of the crossover one would think of 
a scenario in which the correlation length does not play 
a crucial role and we will further comment on this issue 
in the final paragraphs. 

In this respect it is interesting to consider the behavior 
of the susceptibility with respect to external parameters. 
Currently there are great efforts to measure this quantity 
directly [2KH27] the main motivation being the proposed 
existence of a direct connection between xt and the size 
of dynamical heterogeneities Ncorr and thus inferring a 
monotonous increase of the latter. As discussed in |3] 
within SBR it is natural to associate these susceptibil¬ 
ities to d\nTa/da, the resulting plot obtained from the 
3D data is shown in fig. (§ and is similar qualitatively to 
the result of the simplified model. However within SBR 
the increase of dynamical susceptibilities is accompanied 
in the super-cooled region by a decrease of the correla¬ 
tion length and this questions the interpretation of the 
experimental data as evidence of an increase of Ncorr- 
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FIG. 6. Peak of the thermal susceptibility xt as a function 
of relaxation time Tq from SBR theory in 3D (solid line); 
Symbols: experimental data plotted as in Ref. [24] (circles: 
Lennard Jones mixture; squares: hard spheres; BKS silica: 
diamonds; triangles: propylene carbonate; inverted triangles: 
glycerol; open circles: OTP; open squares: salol). Thin line: 
expression proposed in Ref. |24 |. 

From fig. (§ we see that in the supercooled regime 
fluctuations of the B-profile increase constantly in am¬ 
plitude (not in size). We should notice that the theory 
however cannot sustain too large fluctuations and must 
be abandoned at some point. The reason is that the time 
when the total correlator G{t) (controlled by B) enters 


the a-regime becomes much larger than the time when 
the fast regions (corresponding to B[x)/B ^ 1) have en¬ 
tered the a regime. On the other hand as soon as the 
fast regions enter the a regime the description should be 
abandoned because g{x, t) is large (although locally) and 
the scaling condition g{x,t) <C 1 for the validity of the 
critical MCT equation and correspondingly SBR ceases 
to be valid. In particular on time scale where G{t) be¬ 
comes 0(1) and negative the fast regions would have a 
local g{x, t) negative and very large in absolute value, but 
this cannot happens because the density-density correla¬ 
tor F°(q) -b g{x, t)ii{q) must remain positive everywhere. 
What happens when fluctuations become too large can¬ 
not be predicted from SBR. One possibility is that the 
description provided by SBR remains valid because fluc¬ 
tuations are somehow damped by quantitative correc¬ 
tions. Indeed not only the separation parameter but 
also the coupling constants drift with the external pa¬ 
rameters and this changes the amplitude and extension 
of the fluctuations. Another possibility is that one enters 
a full-fledged activated regime and it is tempting to make 
some conjectures on it inspired by SBR. In particular the 
fact that in SBR the correlation length remains relatively 
large, albeit decreasing, below Tc suggest that this ac¬ 
tivated regime could display important differences with 
respect to ordinary activated dynamics, which is driven 
by microscopic {ie. at the single particle scale) events 
occurring exponentially rarely in time but lasting for mi¬ 
croscopic times. In this non-standard activated pictures 
the elementary events (of which the peaks would be the 
precursors) are non-standard in the sense that they in¬ 
volve a relatively large number of particles (although not 
diverging but actually decreasing) and have an intrinsic 
time-scale considerably larger than the microscopic one. 
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